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"j^" ■ Recent progress in exclusive charmless hadronic decays of the B meson is discussed. 

OV 

Q\ , 1. Introduction 

q ! Recently there has been a remarkable progress in the study of exclusive charmless B 

decays, both experimentally and theoretically. On the experimental side, many new two- 
body decay modes were discovered by CLEO M: 

m : 

B^7]'K + , 7]'K° S , ir + K° s , ir°K + , uoK + , uK° s , luit + , (j)K*. (1) 

> 

Moreover, CLEO has improved upper limits for many other channels. Therefore, it is a field 
(N) ■ whose time has finally arrived. On the theoretical aspect, there are two important issues 

to be addressed: (i) the renormalization scheme and scale dependence of hadronic matrix 
elements, and (ii) nonfactorizable effects in charmless B decays. A fascinating progress in 
dealing with the above-mentioned theoretical issues has been made over the last few years. 
In this talk I'll first discuss the theoretical progress and then proceed to elaborate the decay 
B — > r]'K which has received a lot of attention recently. 



2. Renormalization scale and scheme dependence of hadronic matrix element 



The relevant effective weak Hamiltonian for hadronic weak B decay is of the form 

(2) 



x. 



10 

K{cxO\ + c 2 O u 2 ) + A c ( Cl OJ + c 2 O c 2 ) - A t V Ci Oi 



where X u = V u bV* q . The Wilson coefficients Cj(//) in Eq. (2) have been evaluated at the 
renormalization scale fi ~ mj to the next-to-leading order. Beyond the leading logarith- 
mic approximation, they depend on the choice of the renormalization scheme. The mesonic 
matrix elements are customarily evaluated under the factorization hypothesis. In the naive 
factorization approach, the relevant Wilson coefficient functions for color-allowed external 
iy-emission (or so-called "class-I") and color-suppressed (class-II) internal l^-emission am- 
plitudes are given by a\ — C\ + c 2 /N c , a 2 = c 2 + C\/N c , respectively, with N c the number of 
colors. Inspite of its tremendous simplicity, naive factorization encounters two major diffi- 
culties. First, it never works for the decay rate of class-II decay modes, though it usually 
operates for class-I transition. Second, the hadronic matrix element under factorization is 
renormalization scale \i independent as the vector or axial-vector current is partially con- 
served. Consequently, the amplitude C;(/z) (0)f act is not truly physical as the scale dependence 
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of Wilson coefficients does not get compensation from the matrix elements. The first dif- 
ficulty indicates that it is inevitable and mandatory to take into account nonfactorizable 
contributions, especially for class-II decays, to render the color suppression of internal W 
emission ineffective. The second difficulty also should not occur since the matrix elements 
of four-quark operators ought to be evaluated in the same renormalization scheme as that 
for Wilson coefficients and renormalized at the same scale /i. 

To circumvent the aforementioned second problem, one should evaluate perturbative 
QCD and electroweak corrections to the hadronic weak matrix elements parametrized by 
the matrices rh s and rh e , respectively, so that p[] 



Ci(n)(Oi(n)) = Ci(fi) 



I+ ^)^ (/i) + £^ e(/i) 



(Of e ) 



cf (Of e ). 



(3) 



Then factorization is applied to the matrix elements of tree operators; that is, before employ- 
ing factorization to the four-quark operators, it is necessary to absorb all corrections to O tree 
into the effective coefficients c^ ff . One-loop penguin corrections and vertex corrections to the 
operators Oi have been calculated in P, [j[ f|. One can explicitly check that the effective 
Wilson coefficients cf 1 are indeed renormalization scheme independent and approximately 
renormalization scale independent ||. It should be stressed that the effective penguin co- 
efficients cf ff (i — 3, • • • , 10) take into account the effect of penguin diagrams not only with 
top quark exchange characterized by the Wilson coefficients q but also with internal u and 
c quarks induced by the tree operators O" and Of, respectively. 

3. Nonfactorizable effects in charmless hadronic B decay 

As stressed in the last section, it is mandatory to take into account the nonfactorizable 
effects, especially for class-II modes. For B — > PP or VP decays, nonfactorizable contribu- 
tions can be lumped into the effective parameters a± and a 2 pj: 



,cff 



off 



(4) 



where Xi are nonfactorizable terms and receive main contributions from the color-octet cur- 
rent operators. Phenomenological analyses of two-body decay data of D and B mesons 
indicate that while the generalized factorization hypothesis in general works reasonably well, 
the effective parameters a\ 2 do show some variation from channel to channel, especially 
for the weak decays of charmed mesons || 0. An eminent feature emerged from the data 
analysis is that a| ff is negative in charm decay, whereas it becomes positive in bottom decay 



.off' 



(D -> Kn) 0.50 , af(B -> D 



7T 



0.26, 



which in turn implies 

X2(a* ~ m c] D — > Kn) ~ —0.36 , Xvin ~ m b\ B — > Dn) ~ 0.11. 



(5) 



(6) 



The observation |x 2 (-B)| ^ |X2(-D)| is consistent with the intuitive picture that soft gluon 
effects become stronger when the final-state particles move slower, allowing more time for 
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significant final-state interactions after hadronization ||. Phenomenologically, it is often to 
treat the number of colors iV c as a free parameter and fit it to the data. Theoretically, this 
amounts to defining an effective number of colors by 1/Nf = (1/N C ) + \. It is clear from 
Eq. (6) that 

Nf(D -> Kit) > 3, Nf(B -> Du) ~ 2. (7) 



It is natural to ask that does the naive factorization approach also fail in charmless B 
decays ? If so, how large is the nonfactorizable effect ? Since the energy release in charmless 
two-body decays of the B meson is generally slightly larger than that in B — > D^ir, it is 
expected that Nf for the B decay into two light mesons is close to Nf(B — > Dtt) « 2. It 
is pointed out in |H| that the parameters a 2 , a 3 and a 5 are strongly dependent on Nf and 
the rates dominated by these coefficients can have large variation. For example, the decay 
widths of B~ — > uK^~, B° — > uK°, pK*°, B s — > tjuj, T}<f>, uxfi, • - - , etc. have strong iV c 
dependence |IIJ . We have shown recently in [|TTJ that the branching ratio of B~ — > uoK~ has 



its lowest value of order 1 x 10 6 near Nf ~ 3 — 4 and hence the naive factorization with 
Nf = 3 is ruled out by experiment, B{B ± -> cjK ± ) = (l.5i^ ± 0.3) x 10~ 5 M. 



However, it is not easy to discern between Nf = oo and Nf = 2 in B — > uK decays, and 
it becomes important to have a more decisive test on Nf. For this purpose, we shall focus on 
the decay modes dominated by the tree diagrams and sensitive to the interference between 
external and internal H^-emission amplitudes. The fact that Nf = 2 (Nf = oo) implies 
constructive (destructive) interference will enable us to differentiate between them. Good 
examples are the class-Ill modes: B — > tt ^, t]^, ir°p , ■ ■ ■. We found that O] the 



averaged branching ratio of B ± — ► um^ 1 has its lowest value of order 2x10 6 at Nf = oo and 
then increases with l/iV c eff Since experimentally B(B ± -> uit^ = (l.ltoi ± 0.2) x 10~ 5 , 
it is evident that Nf = oo is disfavored by the data. Measurements of interference effects in 
charged B decays B~ — > 7r~(p~)7r°(p ) will help determine the parameter Nf. For example, 
the ratio R x = 2B(B~ -> ir~p )/B(B° -> 7r-p+) is calculated to be 2.50 for Nf = 2 and 0.26 
for Nf = oo. Hence, a measurement of R±, which has the advantage of being independent 
of the Wolfenstein parameters p and 77, will provide a decisive determination of the effective 
number of colors Nf. 

4. Exclusive charmless hadronic B decays to 7/ and rj 

The CLEO collaboration has recently reported the preliminary branching ratios for the 
exclusive decay B — > rj'K dominated by gluonic penguin daigrams [[!]: 

B(B ± -> rfK*) = ~ [B{B + -> ^A^ + ) + -> 7/-FT)] = (7.lj£f ± 0.9) x 10~ 5 , 

-> ?/ K°) = ~ [/3(5° r/if ) + B(B° -> r/A )] = (5.3±|:| ± 1.2) x 10~ 5 . (8) 

Early theoretical estimate of the B ± — > t]'^ 1 * 1 branching ratio [|l|, |1^, || lies in the range of 
(1 — 2) x 10 -5 . [] The CLEO result thus appears to be abnormally large. The question is 

x The prediction B(B ± -> rfK*) = 3.6 x 1(T 5 given in [Q is too large by about a factor of 2 because the 
normalization constant (i.e. l/v3) of the 770 wave function was not taken into account in the form factor 
Fq Vo . This negligence was also made in some recent papers on B — > ry'X. 
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then can the CLEO observation of B — > rj'K be accommodated in the standard model ? Do 
the new data imply new physics ? The theoretical interest and speculation in this subject 
has surged, as evidenced by the recent literature [4,14-20] that offer various interpretations 
on the unexpected large branching ratios. 

In order to illustrate the problem clearly we choose the following parameters for calcula- 
tion: 

F BK (0) = 0.34, VSF o Br,o (0) = 0.254, f = f 8 = f n , m s = 150MeV, 9 = -19.5°, (9) 

and N£ s = 2, where 9 is the 77— rf mixing angle. Using the renormailization scale and scheme 
independent effective Wilson coefficients cf discussed in Sec. II, we find 



B(B ± - r/K*) = l 1Ax 1(r5 ' for v = °- 35 ' p = °- 08 ' 
1 ' ' \l.6xl0" 5 , for 77 = 0.34, p = -0.12. 



(10) 



In the ensuing discussion, we will use (9) and (10) as the benchmarked values to be compared 
with. Since the choice of form factors and light quark masses is uncertain, one may argue 
that the CLEO data (8) can be fitted by choosing a small strange quark mass and/or large 
form factors F BK and Fq Vo [15]. For example, the above model estimate with m s = 55 MeV 



or F BK (0) = 0.63 will fit to the central value of B(B ± — > rfK^). However it is dangerous to 
fit the parameters to a few particular decay modes. The point is that comparison between 
theory and experiment should be carried out using the same set of parameters for all decay 
channels. Indeed the measured branching ratio of B — > nK puts a constraint on the strange 
quark mass and it indicates that m s cannot be too small. In the SU(3) limit we have the 
relation F Bn = F BK . Most of the existing QCD-sum-rule and quark model calculations 
show that F B7r± (0) < 0.33 (for a review, see [p0[]). We shall see below that a severe constraint 
on F Bn (0) can be derived from the current limit on the decay B + — > rj-K + . Since SU(3) 
breaking is expected at most of 30% level, it is very unlikely that Ff K (Q) can deviate much 
from 0.33 . Likewise, the nonet symmetry relation \/3F BrK = F BlT implies that \^3F Br>0 (0) 
cannot be too large than the model estimate, say 0.254 given in (9). In short, the parameters 
given in (9) cannot be modified dramatically without violating SU(3) symmetry relation and 
experimental observation of other decay modes. 

Nevertheless, we can adjust the parameters in (9) slightly to improve the discrepancy 
between theory and experiment. The key point is that an accummulation of several small 
enhancement may eventually lead to a sizable enhancement. First, the current quark mass 
m s = 150 MeV in (9) is defined at the renormalization scale p — 1 GeV. For reason of consis- 
tency, one ought to apply the small running quark mass m s (rrib) — 105 MeV in calculation. 
Second, we use fs/ fn = 1-38 ± 0.22 and fo/ f n = 1.06 ± 0.03 to take into account SU(3) 
breaking effects in decay constants [^2 . Third, for the form factor F B m we follow [|j] to use 



y/3F Br, ° (0) = 0.33, which is slightly larger than the value of 0.254 obtained in Fourth, 



previously we employed the mixing angle 9 = —19.5° so that the wave functions of the r\ and 
7/ have the simple expressions: r\ = -^m{uu + dd — ss) and 7/ = ^(«« + dd + 2ss). Here we 
instead use the value 9 = —22.0° ±3.3° extracted in [p2]] . Applying the new value for each of 
the aforementioned parameters individually, we find that the branching ratio of B ± —>■ rj'K ± 
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is enhanced by 62%, 37%, 19%, and 5%, respectively. Hence, the dominant enhancement 
comes from the running strange quark mass and SU(3) breaking in decay constants. When 
all new parameters are employed, we obtain (see Table I) 

B(B± -> WK*) ~ I 3 ' 87 X 10 ~ 5 ' B(R° -> rt K°) ~ \ 3 ' 7 ° X 10 ~ 5 ' (11) 
^ 7]K ] ~ \ 4.04 x 10- 5 , B[B V K } ~ { 3.67 x 10" 5 , (U) 

for 7] = 0.35, p = 0.08 (upper entry) and r\ = 0.34, p = —0.12 (lower entry), respectively. 
Note that B(B ± — > rj'K^) increases with 1/N* and it becomes as large as order of 6 x 10 -5 
at N% = oo. However, as we have stressed before, it is very unlikely to have a large N* 
in the two-body charmless B decay. Therefore, we see that by adjusting the parameters in 
(9) within some reasonable range, the standard penguin contribution can account for the 
observed decay rate of B° — > rj'K but only marginally for B — > r]'K ± . Nevertheless, the 
current data allow for some new contributions (but not necessarily new physics) unique to 
the rj . Of course, we have to await more new data to sort it out. 

There are several mechanisms which are unique to the rj and may enhance the decay 
rate of B — > rj'K. (i) The b — > sg* penguin followed by the transition g* — > grj via the 
QCD anomaly can in principle contribute to the exclusive decay B — ► g'K. This anomalous 
mechanism was originally advocated to explain the observed large inclusive B — > rj + X 
signal p3|, E4J. However, since this mechanism involves a production of a gluon before 



hadronization, it will not play an essential role in low-multiplicity two-body exclusive decays 
unless the gluon is soft and absorbed in the wave function of the rj . Another possibility is 
that the gluon produced from the penguin diagram and the gluon emitted from the light 
antiquark fuse into the rj [Tj| [TJ|. As the average momentum of the gluon emitted from the 



antiquark is in general less than 1 GeV, it is not clear if perturbative QCD is still applicable 
in this case, (ii) The process b — > s + g*g* — ► s + rj involves two gluon production in the 
penguin-like diagram followed by the ?/-gluon anomalous interaction. The decay b — > s+g*g* 



has been calculated in the literature [25|. It appears that the branching ratio arising from 



this mechanism is less than 1 x 10 ||26|1 . (iii) A new internal jy-emission contribution 
comes from the Cabibbo-allowed process b — > ccs followed by a conversion of the cc pair into 
the rj via two gluon exchanges. This new contribution is important since its mixing angle 
V c bV* s is as large as that of the penguin amplitude and yet its Wilson coefficient a| ff is larger 
than that of penguin operators. The decay constant fjy C \ defined by (Olcy^sclr/) = if^° 'q^ 

has been estimated to be |/^ c ' ) | = (50 — 180) MeV, based on the OPE, large- N c approach 



and QCD low energy theorems 0j. It was claimed in [14J ^7] that \f^\ ~ 140 MeV is 



iff 



needed in order to exhaust the CLEO observation of B — > 7] I K ± and B — > rj + X by the 
mechanism b^cc + s-^g' + s via gluon exchanges. However, a large value of seems 

to be ruled out for three reasons. First, the decay constant f^ u ^ is only of order 50 MeV. 
Second, suppose the pseudoscalar content of cc is dominated by the r\ c . Then from the data 
of J/ip — > r] c 7 and J ftp — > 7/7, one can show that \f^\ > 6 MeV, where the lower bound 
corresponds to the nonrelativitsic quark estimate. (When the relativistic effect of the rj in 
J /if} — > 7/7 is taken into account, \f^\ is larger than 6 MeV.) Even when contributions 



from e.g., r)' c , %, • • • to the cc are included, it is arg ued in J5§ that \f^ c) \ < 40 MeV. Third, 



based on the 777 and 7/7 transition form factor data, the range of allowed is estimated 
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Table I. Averaged branching ratios for charmless B decays to rf and rj, where "Tree" refers to 
branching ratios from tree diagrams only, "Tree+QCD" from tree and QCD penguin diagrams, 
and "Pull" denotes full contributions from tree, QCD and electroweak (EW) penguin diagrams in 
conjunction with contributions from the process cc — > t]q. Predictions are for k 2 = m 2 /2, = 2, 

frf = — 15 MeV, rj = 0.35, p = 0.08 (the first number in parentheses) and r\ = 0.34, p = —0.12 
(the second number in parentheses). 



Decay 


Tree 




Tree+QCD 




Tree4 


QCD+EW 


Full 




Exp. Ill 


R ± — 


ry-J K ± 


1.41 x 10" 


-7 


(a nn 


4 1 81 1 fT 


-5 


(3.87, 


4 041 1 fT 


-5 


(5 48 5 691 10" 


-5 


in i+2.b _i_ q ql i n-5 


R± _ 

/ J 




3.56 x 10" 


-7 


48 


3 ^71 1 rr 


-7 


(3.40, 


^ snl i rr 


-7 


(5.05, 8.68) 10" 


-7 


< 0.8 x 10" 5 


B ± - 


n'K* ± 


2.11 x 10" 


-7 


(2.69, 


1.90) 10" 


-6 


(2.94, 


2.11) 10" 


-6 


(5 90 3 24) 10" 


-7 


< 29 x 10~ 5 


B ± - 




5.25 x 10" 


-7 


(0.92, 


1.53) 10" 


-6 


(1.52, 


2.42)10" 


-6 


(2.49, 3.70) lO- 


-6 


< 24 x 10~ 5 


B ± - 


? 7 ' 7 r ± 


1.94 x 10" 


-6 


(1.13, 


1.06) 10" 


-5 


(1.12, 


1.05) 10" 


-5 


^l. 29, 1.21) 10- 


-5 


< 4.5 x 10~ 5 


B ± - 




4.93 x 10" 


-6 


(9.57, 


6.02) 10" 


-6 


(9.82, 


6.24) 10" 


-6 


(1.04, 0.67) 10- 


-5 


< 0.8 x 10" 5 


B ± - 


v'p ± 


3.95 x 10" 


-6 


(1.08, 


1.84)10" 


-5 


(1.08, 


1.84)10" 


-5 


(1.01, 1.71) lO- 


-5 




B ± - 


i 1 p ± 


9.72 x 10" 


-6 


(1.21, 


1.71)10" 


-5 


(1.19, 


1.66)10" 


-5 


ll. 19, 1.63) 10" 


-5 




B d ^ 




5.09 x 10" 


-9 


(3.82, 


3.80) 10" 


-5 


(3.70, 


3.67) 10" 


-5 


(5.22, 5.19) 10" 


-5 


(5.3+2^±1.2)10- 5 


B d ^ 




1.93 x 10" 


-8 


(1.23, 


0.77) 10" 


-7 


(1.97, 


3.11) 10" 


-8 


(2.49, 3.28) 10" 


-7 




B d ^ 


V 'K*° 


3.88 x 10" 


-9 


(2.23, 


2.08) 10" 


-6 


(2.47, 


2.32) 10" 


-6 


(2.44, 2.96) 10- 


-7 


< 4.2 x 10~ 5 


B d ^ 


rjK*° 


1.62 x 10" 


-8 


(5.71, 


6.64) 10" 


-7 


(1.23, 


1.36) 1Q- 


-6 


(2.26, 2.44) 10- 


-6 


< 3.3 x 10~ 5 


B d ^ 


V 'tt° 


3.31 x 10" 


11 


(4.83, 


6.80) 10- 


-6 


(4.63, 


6.52) 10" 


-6 


(5.32, 7.35) 10- 


-6 


< 2.2 x 10" 5 


B d ^ 


r)TT° 


6.70 x 10" 


-9 


(2.71, 


3.67) 10" 


-6 


(2.66, 


3.60) 1Q- 


-6 


(2.86, 3.83) lO- 


-6 




B d ^ 


V'P° 


1.21 x 10" 


-7 


(2.07, 


3.33) 10" 


-6 


(2.01, 


3.24) 10- 


-6 


il. 79, 2.93) lO- 


-6 




B d ^ 


r)p° 


3.46 x 10" 


-7 


(1.28, 


2.20) 10" 


-6 


(1.16, 


1.98) 10- 


-6 


ll. 11, 1.89) 10" 


-6 


< 8.4 x 10" 5 



to be -65 MeV < ffi < 15 MeV |9 



From Table I it is clear that for = —15 MeV, which is consistent with above- 
mentioned constraints, the agreement between theory and experiment for B — > rj'K is 
substantially improved in the presence of large charm content in the rj'. We conclude 
that no new physics is needed to account for the CLEO data of B — > rj'K. We have 
also calculated the branching ratios of other exclusive charmless B decays involving rj' 
and r] (see Table I), where use of = — tanfl/^/ has been made. Three com- 
ments are in order, (i) The effect of cc conversion into the rj' contributes destructively 
to B — > rj'K*. Consequently, the branching ratio of B — > r/'K* is suppressed Q and 
B(B -> r]'K*)/B(B -> rjK*) ~ C^ICT 1 ). If B -> rj'K is assumed to be entirely accom- 
modated by any of aforementioned new mechamisms, the decay rate of B — > rj'K* will 
be predicted to be the same order of magnitude as B — > rj'K |18j. (ii) Contrary to 
B -> 7]{r]')K* decays, we see from Table I that B(B -> rfK)/B(B -> r]K) ~ C(10 2 ) due 
to the destructive interference in the penguin diagrams of B — > r]K. (iii) The electroweak 
penguin effects are in general very small, but they become important for B — > r]K and 



Jcc) 



! ln the two mixing angle parametrization scheme given in p{J], the decay constant fjf^ is much smaller: 



tan 6»i f { J c> with 6»i = -(6° - 9°). 



Tn 



we have employed f^, 



(cc) 



-50 MeV. In that case, B — > rj'K* is dominated by the process cc — > r/o 



and its branching ratio is of order 10 
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B — > r]K* decays due to a large cancellation of QCD penguin contributions in these decay 
modes. 

For B — > i]'(r])7r(p) decays, the mechanism of cc — > 770 is much less dramatic since it does 
not gain mixing-angle enhancement as in the case of B — > r]'(rj)K(K*). Their branching 
ratios are sensitive to the light quark masses m u , and form factors such as F^ n [||]. 
The current experimental limit on the decay £P — > rjir^ puts useful constraints on m q and 
Fq*. The predicted values presented in Table I are for m u {m h ) = 5 MeV, m d {m b ) = 10 
MeV and _F S7r± (0) = 0.30. We find that even a slight increase of F S7r (0) or decrease of m q , 
say m u (mj)) ~ md{mb)/2 ~ 3 MeV, will make the decay rate of 5^ —» rpr 1 * 1 exceeding the 
present upper bound significantly. We also see that a negative p, which in turn implies a 
unitarity triangle 7 in the range 90° < 7 < 180°, is preferred [ H| . By contrast, the present 



experimental value of the ratio R2 = T(B° — ► 7r =F /i' ± )/r(i? ± — > 7r ± A"°) favors a positive p 
Q. Note that a positive p is also preferred by the limit on the ratio AM S /AM^ |30 |. Clearly 
more data of B ± — > r/^ and i?2 are needed to pin down the sign of p. 
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